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ABSTRACT 

The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point 
Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 
10° < |Z| < 65° and |6| < 1°, was used to determine the distribution of 
stars in Galactic longitude, /, latitude, b, and apparent magnitude, m. The 
counts versus longitude can be approximated by the modified Bessel function 
= Nq{1/Iq)Ki{1/Iq), where Iq is insensitive to limiting magnitude, band choice, 
and side of Galactic center: Iq = 17° — 30° with a best fit value in the the 4.5/im 
band of Iq = 24° ± 4°. Modeling the source distribution as an exponential disk 
yields a radial scale length of = 3.9 ± 0.6 kpc. 
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There is a pronounced north-south asymmetry in source counts for |Z|^30°, 
with ~25% more stars in the north. For / = 10 — 30°, there is a strong enhance- 
ment of stars of m = 11.5 — 13.5 mag. A hnear bar passing through the Galactic 
center with half-length Rbar = 4.4 ± 0.5 kpc, tihed by </> = 44° ± 10° to the 
Sun-Galactic Center line, provides the simplest interpretation of this data. We 
examine the possibility that enhanced source counts at / = 26 — 28°, 31.5 — 34°, 
and 306 — 309° are related to Galactic spiral structure. Total source counts are 
depressed in regions where the counts of red objects {rriK — ?«[8.o] > 3) peak. In 
these areas, the counts are reduced by extinction due to molecular gas and/or 
high diffuse backgrounds associated with star formation. 

Subject headings: infrared: general - infrared:ISM - infrared:stars - survey - 
stars:general - ISM:general - Galaxy:structure - Galaxy:stellar content 

1. GLIMPSE and the stellar structure of the Galaxy 

Since the pioneering work of Kapteyn (1922), star counts have been an important avenue 
in exploring the stellar structure of the Galaxy. However, progress on understanding the 
stellar content of the inner Galactic disk has been impeded by dust obscuration of visible 
light from stars. With the advent of infrared detectors and surveys of increasingly high 
sensitivity and resolution, significant progress has been made, leading, for example, to strong 
evidence for a Galactic bar (reviewed by Gerhard 2002). 

GLIMPSE , a Spitzer Legacy Science Program, has imaged the Galactic plane from 
Galactic longitude^ |/| = 10° to 65° and Galactic latitude \b\ < 1° at 3.6, 4.5, 5.8, and 
8.0 /xm using the IRAC (Infrared Array Camera, Fazio et al. 2004) on the Spitzer Space 
Telescope (Werner et al. 2004). Angular resolutions range from 1.4" at 3.6 /xm to 1.9" at 
8 //m. We present the mid-infrared distribution of stars in the Galaxy as determined by 
GLIMPSE . Properties of the GLIMPSE Point Source Catalog and Archive are summarized 
in Table 1. Observing strategy, details of source selection, quality flags, and photometric 
accuracy are provided in the GLIMPSE Science Data Products document (Meade et al. 
2005), the GLIMPSE Quahty Assurance Document (Churchwell et al. 2005), and Benjamin 
et al. (2003). 

We have compared the results obtained for all four /i?>lC bands, explored the effect 
of using the Archive instead of the Catalog, explored the sensitivity of our results to our 



|Z| < 180° is the angular distance from Galactic center. 



-3- 



faint magnitude limit, divided the survey region into three latitude strips of Ab — 0.3° and 
considered each strip separately, and experimented with the effect of removing sources with 
various quality flags. The parameters of all the principal features discussed below remained 
unchanged. 

2. GLIMPSE Counts as a Function of Longitude and Latitude 

The Catalog and Archive sources were binned into data cubes of 0.1 "(longitude) x 
0.1° (latitude) x 0.1 magnitude for all four IRAC hands. The average number of sources in 
each 0.1°x 0.1° bin is ~1,400 for the Catalog and ~2,200 for the Archive. Figure 1 shows 
the average number of 4.5;Lim sources per square degree in the magnitude range m = 6.5 
to m = 12.5. The curves indicate the source-count weighted Galactic latitude distribution, 
bcen{l) = ^^N{n) ' '^hich ranges from been = —0.05° to +0.05°. As |/| increases from 10° 
to 65°, the counts decrease by a factor of ~ 12, but there is a pronounced asymmetry 
between the north and the south in the inner Galaxy. From 10°-22° there are ~ 25% 
more sources than in the equivalent southern longitude range. There are also apparent 
enhancements in the counts for Z = 26° - 28°, 31.5° - 34°, and 306° - 309°. Figure 1 also 
shows contours of the surface density of ~ 40, 000 unusually red {uik — ""^[s.o] > 3.0) sources 
using the GLIMPSE and 2MASS Catalogs (Cutri et al. 2001). These objects are strongly 
anti-correlated with the total GLIMPSE counts. The regions of depressed counts coincide 
with high column densities of ^^CO (Simon et al. 2001; Shah 2005, priv. comm.) or thermal 
radio continuum emission (i.e. high mass star formation). 

Figure 2 shows that the number of sources per square degree brighter than magnitude 
m is well-fit by the first order modified Bessel function of the second kind, A^(band, m, |/|) = 
A'o(|/|//o)-^i(KI/^o), which is expected for an exponential disk population (Merrifield, priv. 
comm.). Our "best-fit" curve excludes data interior to \l\ — 35°, from I — 306.5 — 309° 
(enhancement) and I — 301.5 — 306.5° (missing data). The inner boundary was varied 
from \l\ — 25° to \l\ — 45° to estimate systematic uncertainties due to this choice. The 
range for all the IRAC bands is Iq — 11 — 30°. The best fit in the 4.5/im band is Iq — 
24.2 ± 0.3°(random) ± 3.6° {2a systematic); the slope and level of the fits for the north 
and south are in marginal agreement. Comparison of this fit to numerical models of a disk 
population of M and K giants, which are assumed to dominate the source counts (Wainscoat 
et al. 1992), yields an exponential scale length of iJ* = 3.9 ± 0.6 kpc. 

For longitudes greater than |/| = 30° and excluding the enhancements mentioned above, 
the fitting functions for the 4.5/um band predict the observed source densities to within 20% 
for a faint magnitude cutoff of 8 and 10 mag, and to within 10% for faint magnitude cutoffs 
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of 12 and 14 magnitudes. Interior to \l\ — 22°, the fits overpredict the observed counts by as 
much as 30% in both the northern and southern plane. This is not due to confusion as the 
drop is independent of the magnitude cut; even at faint magnitudes, we are not affected by 
confusion (Churchwell et al. 2005). 



3. GLIMPSE Counts as a Function of Magnitude 

Figure 3 shows the number of sources as a function of magnitude for three north/south 
pairs of directions. Figures 3 and 4 were used to determine the "effective" sensitivity hmit, 
nisens) as a function of longitude given in Table 1. The curves shown in Figure 3 can be fit by 
\ogiQn{m) = am + b. Since m = — 2.5 logiQ(S'/S'o), where S is the flux density at magnitude 
m and 5*0 is the flux density for zero magnitude, we can write N{S) = Nq (^-§^^ , where 
iVo = (2.5/5*0 In 10)10^ sources deg^^ Jy^^ and a = 2.5a + 1 . The average a over the 
magnitude range 8 < m < 11 lies in the range 1.95 to 1.83 for the entire survey area. 

Figure 3 shows a "hump" at magnitudes fainter than m = 12 in the inner Galaxy 

towards the north but not the south. This hump is seen in both the Catalog and Archive 
data in the 3.6yum and 4.5/xm bands. It is not detected in the 5.8;um and 8.0/im bands 
due to lower sensitivity. To study this further, in Figure 4 we plotted the power-law slope, 
a = 2.5^^^^^ + 1 (with Am = 0.1) of the source distribution as a function of Galactic 
longitude and magnitude. The hump persists over the range / =10-22° and / =24-29°. It 
seems hkely that the gap at I ~ 23° is caused by high extinction or high diffuse flux. The 
central magnitude of the hump at I — 10° is nih — 12.5 ± 0.1 mag with a dispersion of 
afi — 0.25 ± 0.10 mag; the amplitude at the peak is 20% above the power- law fit. The 
shift of the central magnitude of the hump to brighter magnitudes with increasing Galactic 
longitude is given by druh/dl — —0.025 ± 0.005 mag deg~^. 



4. Interpretation and Comparison with Previous Results 

We now summarize the principal features we have found and possible implications for 
our understanding of Galactic structure. We assume the distance to the Galactic center is 
Ro — 8.5 kpc. 

Exponential Disk: For |Z|^22°, the counts in all four IRAC bands can be fit by A'"(band, m, |Z|) = 
-^o(|^|/^o)-f^i(|^|/^o)- For the 4.5 jjm band, the best fit (excluding regions of obvious enhance- 
ment) yields Iq — 24.2±3.6°(systematic). For all four IRAC bands and for several magnitude 
cutoffs, we find that Iq — 17° — 30°. This can be compared with our results obtained by 
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binning the 2MASS K band sources in an identical fashion which yields Igc > 40°. Presum- 
ably this larger number is due to higher extinction in the K band. Converting this angular 
scale length to a radial scale length requires modeling of the spatial distribution of stars, the 
infrared luminosity function, and extinction (Cohen 1993). Preliminary modeling indicates 
that the majority of our sources are M giants and late K giants, and yields a radial scale 
length for the exponential disk of if* = 3.9 ± 0.6 kpc. Interior to |/| = 22° there are ~ 
20% fewer sources than predicted by this fit, which may indicate a stellar "hole" (or more 
accurately, deficit) interior to Galactocentric radius, Rh — 3.2 ±0.3 kpc. This is comparable 
to Rh = 2.9 - 3.3 kpc derived from COBE/DIRBE data (Preudenreich 1998). 

Galactic Bar: We see a strong north/south asymmetry in counts for |/|^30°, with ~25% 
more sources in the north and a distinct hump in the number of stars at apparent magnitude 
m ~ 12.5 over the longitude range I — 10 — 30°. Both of these facts suggest that we have 
detected the Galactic bar. Estimates for the bar orientation (relative to the line connecting 
the Sun to the Galactic center) typically range from (p = 15° — 35° (Gerhard 2002) . If the 
"hump" seen in Figure 3 and 4 is due to stars in the bar with a well-defined peak at constant 
absolute magnitude, we can estimate the bar orientation and length. In Figure 4, we have 
plotted the magnitude and longitude expected for a linear bar in the plane of the Galaxy 
passing through the Galactic center with a foreground extinction a[4,5]((i) = 0.05 mag kpc~^, 
using the fiducial av{d) = 1 mag kpc"^(Mihalas & Binney 1981), Av/Ak ~ 8.8 — 7.5 
(CardeUi, Clayton, & Mathis 1989) and ^[4.5]/^^ = 0.43 (Indebetouw et al. 2005). The 
best fit to the data yields = 44 ± 10°, R^ar = 4.4 ± 0.5 kpc, and M[4.5] = -2.15 ± 0.2 
mag; the effect of increasing extinction is to increase and M. Pitting the data with zero 
extinction yields M[4.5] = —1.8 mag. The requirement that M[4.5] = —2.15 mag to match 
the data suggests that the hump is due to K2-3 giants (Cohen 1993; Hammersley et al. 
2000). Figure 4 shows that the predicted nihump ioi the model bar falls below our sensitivity 
limit at / = 350°, but the observed increase in counts seen in Figure 1 from / =345° to 
/ =350° is consistent with our derived bar length. Although a bar seems to be the simplest 
interpretation of these results, other configurations should be explored. 

Stellar Enhancements/Spiral Arms(?): Outside / = 22°, we find three regions where 
counts are 20% higher than the exponential fitting functions. These are at / = 26 — 28°, 
31.5 — 34°, and 306 — 309°. It is tempting to associate these features with spiral arm 
tangencies. The presence of the "hump" in the I — 26 — 29.5° region argues that it is a 
continuation of the same structure observed from I — 10 — 20°, complicated by the presence 
of extinction. Unresolved questions regarding extinction also affect interpretation of the 
31 — 34° feature. The third feature could be identified with a previously known spiral arm 
tangency (Centaurus, / ~ 309°). No plausible enhancement is star counts is seen towards 
the Sagittarius spiral arm tangency at I ~ 49°. 
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Extinction/ Diffuse Background: There is a clear anticorrelation between the GLIMPSE 
Catalog counts and the counts of sources with rriK — m\^?,.o\ > 3.0 in certain regions (shown 
in Figure 1). Counts can be depressed by cither increased extinction or high diffuse back- 
ground, both of which are associated with star forming regions and molecular clouds. It 
is clear that the effects of extinction and diffuse emission on the source counts deserves fu- 
ture scrutiny. For this reason, the parameters of the structures discussed here should be 
considered preliminary. 

Many of the principal structures seen in this work have been detected or hinted at 
in previous infrared investigations. What differentiates the GLIMPSE results from these 
previous works is the uniform and complete sampling of the Galactic plane at ^2 arsec 
resolution, the significantly reduced extinction in the mid-infrared, and the very large sample 
of GLIMPSE sources. Even without any color-selection of our sources, we find that we are 
able to detect and confirm several fundamental Galactic structures. 

We thank the referee, Stephan Jansen, Mike Merrifield, Paul Rybski, and Ronak Shah 
for help and suggestions, and acknowledge support for this program through NASA contracts 
1224653 (UW-W), 1224653 (UW-M), 1225025 (BU), 1224681 (UMd), 1224988 (SSI), 1259516 

(UCB), 1253153 (UMn), 11253604 (UWy) by the Jet Propulsion Laboratory. Caltcch under 
NASA contract 1407. We also gratefully acknowledge use of data products from the Two 
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Fig. 1. — Logio(sources deg~^) of 4.5 /iin band sources from the GLIMPSE Point Source Cat- 
alog in the magnitude range 6.5 < m < 12.5 as a function of Galactic latitude and longitude. 
Sources have been binned by 0.1° (longitude) x 0.1°(latitude). Notable features include 
a decline in the number of sources with distance from Galactic center and a pronounced 
north-south asymmetry in the inner Galaxy between the north(/ = 10 — 22°) and south 
(Z = 350 — 338°). The number of sources is nearly symmetric about the Galactic midplane; 
the solid line shows the position of the source-count weighted average latitude. The two 
contours, logio(sources deg~^) = 2.8, 3.3 show the surface density of red {rriK — ?^[8.o] > 3) 
sources; note the anti-correlation with the total counts. These directions are labeled with 
corresponding star forming regions or R for strong radio continuum (Altenhoff et al. 1970; 
Haynes, Caswell, & Simons 1978). The dark stripe from / = 301 — 306° is a region of missing 
data. 
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Fig. 2. — Average number of sources per square degree detected by GLIMPSE as a function 
of longitude in both the northern Galactic plane (red curves, coordinates on lower axis) and 
southern Galactic plane (blue curves, coordinates on upper axis). Average source densities 
were obtained by binning the stars into strips of 0.1 degree in longitude and 1.8 degrees in 
latitude. Fits to both curves (see text) over the range I = 35° — 65° (excluding regions of 
clear enhancements) are also shown. The four sets of curves show the effect of changing 
the faint mangitude cutoff. For each drop in magnitude, the number of sources increases 
by a factor of 3.5 (at 9 mag) to 2.3 (at 13 mag). The vertical lines indicate the estimated 
directions of spiral arms using different tracers (Englmaier & Gerhard 1999). The solid line 
shows their adopted values; the dotted lines indicate the range of published estimates. 
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Fig. 3. — Number of sources from the GLIMPSE Point Source Catalog as a function of mag- 
nitude for three northern and three southern directions. The outer galaxy (/ = 55.5°// = 
304.5°) and middle galaxy (/ = 35.5°// = 324.5°) curves have approximately the same am- 
plitude and slopes. The inner galaxy (/ = 15.5°// = 324.5°) shows a significant north/south 
asymmetry; the northern direction also shows a bump in source counts at a magnitude of 
m ~ 12.2. For all directions shown, the number of sources has been averaged over a 1.0 
(longitude) x 1.8 (latitude) region. The Catalog is truncated at 6.5 mag due to detector 
nonlinear ity for sources brighter than this limit. 
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Fig. 4. — Power-law exponent of counts as a function of flux density plotted as a function 
of apparent magnitude and Galactic longitude. The position of the ~ 12 mag hump seen 
towards I = 15.5° in Figure 3, is seen here to vary consistently in both longitude and 
magnitude. The locus of magnitude and longitude of a model bar consisting of stars of 
absolute magnitude M[4.5] = —2.15, foreground extinction a[4.5](r) = 0.05 mag kpc~^ and 
three different position angles 0, are shown in black in the upper panel and in white in the 
lower panel. The dots indicate R =3, 4, and 5 kpc points along the bar. The dotted hues 
show the same position angles for zero extinction, with M[4.5] = —1.8. 



- 12 - 



Table 1. GLIMPSE Catalog/ Archive Source Information 



IRAC 
Band 


(/xm) 


(Jy) 


Ak 


rrisef 
mag 


mag 


mag 


Catalog Sources^ 
North South^ 


Archive Sources' 
North South 


1 


3.55 


277.5 


0.56 ±0.06 


14.2 


7.0 


13.3-13.6 


14.775 


14.255 


21.420 


22.044 


2 


4.49 


179.5 


0.43 ±0.08 


14.1 


6.5 


13.3-13.6 


14.768 


14.250 


19.797 


19.423 


3 


5.66 


116.5 


0.43 ± 0.10 


11.9 


4.0 


11.7-12.3 


5.768 


5.291 


6.095 


5.594 


4 


7.84 


63.13 


0.43 ± 0.10 


9.5 


4.0 


11.0-12.4 


4.426 


3.959 


4.749 


4.268 



^Vega isophotal wavelengths and IRAC zero magnitudes from Cohen (2005, priv. communication) 
^Extinction from Indebetouw et al. (2005) 

^ GLIMPSE Point Source Catalog selection limits and brightness cutoff limits from Meade et al. 
(2005) 

"^Millions of sources 

"Effective" Catalog sensitivity limit varies over the longitude range |/| = 10° to |/| = 65° 
^Southern Catalog/ Archive is still missing ~1% of survey area 



